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We investigate how quantum correlations (quantum discord (QD)) of a two-qubit one dimensional 
XYZ Heisenberg chain in thermal equilibrium depend on the temperature T of the bath and also on 
an external magnetic field B. We show that the behavior of thermal QD differs in many unexpected 
ways from thermal entanglement. For example, we show situations where QD increases with T when 
entanglement decreases, cases where QD increases with T even in regions with zero entanglement, 
and that QD signals a quantum phase transition even at finite T. We also show that by properly 
tuning B or the interaction between the qubits we get non-zero QD for any T and we present a new 
effect not seen for entanglement, the "regrowth" of thermal QD. 
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Introduction. Since the seminal work of John S. Bell 
PQ, who brought to the realm of experimental Physics 
the ideas of Einstein, Podolsky, and Rosen [2 concern- 
ing the non-local aspects of Quantum Mechanics, it be- 
came clear that the constituents of some quantum com- 
posite systems possessed correlations among themselves 
unachievable in the classical world. Those states belong 
to the class of entangled states although not all entangled 
states possess stronger-than-classical correlations in the 
sense of violating Bell inequalities [3 j. However, given a 
composite quantum state whose constituents are corre- 
lated, how can one tell the origin of the correlations? In 
other words, how can one divide the total correlation into 
a classical part and a purely quantum one? This is par- 
ticularly important for mixed states since their quantum 
correlations are many times hidden by their classical cor- 
relations. An answer to these questions is given by the 
quantum discord (QD), a measure of the quantumness 
of correlations introduced in Ref. [4 . QD is built on 
the fact that two classical equivalent ways of defining 
the mutual information turn out to be inequivalent in 
the quantum domain. In addition to its conceptual role, 
some recent results [5] suggest that QD and not entan- 
glement may be responsible for the efficiency of a mixed 
state based quantum computer. 

Due to its fundamental and practical significance we 
wish to investigate the amount of QD in a concrete sys- 
tem such as a pair of qubits (spin- 1/2) within a solid 
at finite temperature, whose interaction is given by the 
Heisenberg model (XYZ model in general). Such Heisen- 
berg models can describe fairly well the magnetic prop- 
erties of real solids [6 and is well adapted to the study 
of the interplay of disorder and entanglement as well 
as of entanglement and quantum phase transitions [7]. 
We characterize the dependence of QD on the tempera- 
ture and also on external magnetic fields applied to the 
qubits. We also compare the behavior of QD against 
the entanglement of formation (EoF) between the two 
qubits [8HTU]. We obtain several interesting and novel 
results for the behavior of QD at non-zero T, many of 



them in contrast to the behavior of EoF. First, we show 
that QD can increase with temperature in the absence 
of external fields applied to the qubits. This is in sharp 
contrast to the behavior of EoF since one can show [10] 
that this effect never occurs without the presence of an 
external magnetic field. We also show that there exist re- 
gions where QD is different from zero while EoF is always 
zero, confirming a generic feature of QD [11 . In partic- 
ular, we show that for the isotropic XXX model, both 
the ferromagnetic and anti-ferromagnetic Hamiltonians 
possess relatively high values of QD while EoF is absent 
for the ferromagnetic model [8]. These properties of QD 
can have important practical consequences for the real- 
ization of a quantum computer at finite temperatures. 
Indeed, once it is established that what is at stake for 
the correct functioning of a quantum computer is the ex- 
istence of a certain level of QD (just having non-null QD 
might not be enough since almost all states have QD> 
[TT]), we have shown that Heisenberg solids are better 
than previously thought as good candidates for the con- 
struction of a quantum computer at T > 0: for some 
parameter sets they have high values of QD while having 
no EoF at all. We also show that by properly adjusting 
the coupling constants [9] one can achieve any desired 
level of QD for any T. Finally, we show that for finite T 
we observe the sudden-change of QD [12 as we change 
the coupling constant of the Hamiltonian and we show 
a new effect, namely the "regrowth" of the thermal QD 
with increasing T. 

The thermalized Heisenberg system. The Hamiltonian 
of the XYZ model with an external magnetic field acting 
on both qubits is 

H = B(S\ + S z ) + J x SlSl + JySySy + J Z S\S 2 Z , (1) 

with J Xl J y , and J z the coupling constants, S 3 xy z = 
a x,y,z/2, °~i->°~yi an d o~ z the usual Pauli matrices acting 
on qubit j, and B the external magnetic field. We have 
assumed H = 1. The density matrix describing a system 
in equilibrium with a thermal reservoir at temperature 
T (canonical ensemble) is p = exp(—H/kT)/Z, where 
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Z = Tr {exp (— H/kT)} is the partition function and k 
is Boltzmann's constant. Therefore, Eq. leads to the 
following thermal state in the standard basis 



1 

Z 



/An A 12 \ 

Bu B 12 

B 12 B u 

\ Ai 2 A 22 J 



(2) 



Here A u = e~ a (cosh(/3) - 45 sinh(/3)/^), A i2 = - A 
e~ a sinh(/3)/r7, A 22 = e~ a (cosh(/3) +45 sinh(/3)/r/), 
5n = e a cosh(7), Bi 2 = — e a sinh(7), and Z = 
2 (exp (—a) cosh(/3) + exp (a) cosh(7)), where A = J x — 
J y , Z = J x + J y , T] = VA 2 + 165 2 , a = J*/(4fcT), 
/? = 77/(4jfeT), and 7 = £/(4fcT). 

Entanglement. For a pair of qubits there exists an ana- 
lytical expression to quantify its amount of entanglement 
called Entanglement of Formation (EoF) [13] . Given the 
density matrix p describing thermalized two qubits, EoF 
is the average entanglement of the pure state decompo- 
sition of p, minimized over all possible decompositions, 
EoF(p) = mmJ2kPk^(M, where J2kPk = 1> < Pk < 
1, and p = ^2 k Pk \<t>k) (<t>k\- E(<fik) is the entanglement of 
the pure state \<pk) [H]- For a pair of qubits Wootters [15] 
has shown that EoF is a monotonically increasing func- 
tion of the concurrence C (an entanglement monotone), 
EoF = -/(C)log 2 /(C) - (1 - /(C))log 2 (l - /(C)), 
where /(C) = (1 + Vl - C 2 )/2. The concurrence is 
simply [T5] C = max{0, Ai— X 2 — A3 — A4}, where Ai, A 2 , A3, 
and A4 are the square roots of the eigenvalues, in de- 
creasing order, of the matrix R = pp. Here p is the time 
reversed matrix p = (cr* (g) <r 2 ) p* (<j^ (g) a 2 ) . The sym- 
bol means complex conjugation of the matrix p in the 
standard basis {|00) , |01) , 1 10) , |H)}. For a density ma- 
trix in the X-form above C = 2max{0, Ai, A 2 }/Z, with 
Ai = |^i 2 | - VAnA 22 and A 2 = \A 12 \ - B tl . 

Quantum Discord. In classical information theory 
(CIT) the total correlation between two systems (two 
sets of random variables) A and B described by a joint 
distribution probability p(A, B) is given by the mutual 
information (MI), 



I(A, B) = H(A) + H(B) - H(A, B), 



(3) 



with the Shannon entropy H(-) = — Pjlog 2 Pj. Here 
Pj represents the probability of an event j associated to 
systems A, 6, or to the joint system AB. Using Bayes's 
rule we may write MI as 



1(A,B)=H(A)-H(A\B), 



(4) 



where H(A\B) is the classical conditional entropy. In 
CIT these two expressions are equivalent but in the quan- 
tum domain this is no longer true [4 . The first quan- 
tum extension of MI, denoted by 1>(p\ is obtained di- 
rectly replacing the Shannon entropy in Q with the von 
Neumann entropy, S (p) = — Tr(plog 2 p), with p, a den- 



sity matrix, replacing probability distributions. To ob- 
tain a quantum version of Q it is necessary to gener- 
alize the classical conditional entropy. This is done rec- 
ognizing H(A\B) as a measure of our ignorance about 
system A after we make a set of measurements on B. 
When B is a quantum system the choice of measure- 
ments determines the amount of information we can ex- 
tract from it. We restrict ourselves to von Neumann 
measurements on B described by a complete set of or- 
thogonal projectors, {n^}, corresponding to outcomes j. 
After a measurement, the quantum state p changes to 
Pj ~ [(" ® Hj) p (0 (g) Ilj)] /pj, with D the identity opera- 
tor for system A and pj = Tr[(D <g> H J -)p(D <g> IT,-)]. Thus, 
one defines the quantum analog of the conditional en- 
tropy as S (p \ {n^}) = pjS (pj) and, consequently, 
the second quantum extension of the classical MI as 
HI J(P I { n j}) = S (p A ) - S(p I {n^}). The value of 
J (p I {n^}) depends on the choice of {Ilj}. Hender- 
son and Vedral [16] have shown that the maximum of 
J (p I {n^}) with respect to {Hj} can be interpreted as 
a measure of classical correlations. Therefore, the differ- 
ence between the total correlations X (p) and the classical 
correlations Q(p) = sup{ n j. J (p | {n^}) is defined as 



D(p)=l{p)-Q{p), 



(5) 



giving, finally, a measure of quantum correlations [4] 
called quantum discord (QD). For pure states QD re- 
duces to entropy of entanglement [14 , highlighting that 
in this case all correlations come from entanglement. 
However, it is possible to find separable (not-entangled) 
mixed states with nonzero QD [4 , meaning that entan- 
glement does not cause all nonclassical correlations con- 
tained in a composite quantum system. Also, QD can 
be operationally seen as the difference of work that can 
be extracted from a heat bath using a bipartite system 
acting either globally or only locally [17]. For more in- 
sights into QD the reader is referred to [18 . In this work, 
when B — the density operator Q is such that QD is 
given by an analytical expression obtained in [19] . When 
5/0 we computed QD numerically [20] . 

Results. Let us start presenting the important result 
that QD increases with temperature without an exter- 
nal field acting on the qubits (B = 0). This effect 
can be clearly seen when we deal with the XXZ model 
(J x = J y = J and J z 7^ 0). Looking at Fig. [I] panels a 
and 6, we see that this effect happens for several config- 
urations of coupling constants being, thus, dense around 
this region. We should note that such behavior can be 
found for other models than the XXZ as well. As we men- 
tioned before, in the absence of external fields [10] such 
increase with T only occurs for QD and not for EoF. Note 
also that QD starts at zero and then increases with T. 
In particular, for the set of coupling constants shown in 
Fig. [I] the EoF is always zero [10 . Furthermore, one can 
show that the classical correlation decreases when QD 
increases (See Fig. [1]^). We have, therefore, a genuine 
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FIG. 1: (Color online) (a) and (b) QD as a function 
of the absolute temperature kT for the XXZ model with 
B = 0. (a) Here J z — —0.5 and from bottom to top J — 
0.1,0.2,0.3,0.4; (b) Now we fix J = 0.4 and from bottom to 
top J z = -0.8,-0.7,-0.6,-0.5. (c) QD (solid/red line) and 
classical correlations, Q(p), (dashed/black line), as functions 
of kT for J = 0.4 and J z = —0.5. Here and in the following 
graphics the quantities plotted are dimensionless. 



increase of quantum correlations with decreasing classi- 
cal correlation. This non trivial and unexpected effect 
indicates the robustness of quantum correlations against 
classical correlations as we increase the temperature. 

We now move to the XXX model (J x = J y = J z = J) 
with no field. Looking at Fig. [2] we see that although 
EoF is zero in the ferromagnetic region (J < 0) we have 
non-zero values for QD. Also, both in the ferro and an- 
tiferromagnetic regions QD increases if we increase the 
absolute value of J. Actually, for any finite T we can 
find a J big enough such that QD ^ 0. This point is 
justified observing that for any finite value of T, the den- 
sity operator Q when J — >> oo is given by the Bell state 
p = |^)(^| with |^) = \ (|01) - |10>). Besides, the 
density operator in the limit J —> — oo for a finite T is 
a mixed state p = \ (|00) (00| + |11) (11| + \<j>) {(j>\) with 
= i (|01) + 1 10>). The value of QD in this case is 1/3 
and, as expected, the EoF is zero. We should also men- 
tion that QD behaves qualitatively differently for T > 
whether J < or J > 0, achieving the value zero at 
the critical point J = 0. This suggests that QD can sig- 
nal the quantum phase transition (QPT) taking place at 
T = (when we tune J driving the system from its anti 
to its ferromagnetic phase) even at finite T. Note that 
EoF does not signal this QPT at finite T: EoF goes to 
zero before the critical point for T > while QD = 
exactly at the critical point, and only there. However, 
just as was done for T = [21], a detailed study of the 
behavior of QD and QPT for T > is needed. 

Focusing now on the XYZ model we investigate how 



FIG. 2: (Color online) EoF (a) and QD (b) as a function of 
the absolute temperature kT and the coupling J. Both plots 
for the XXX model with B = 0. 

QD behaves while we change the coupling constants for 
a fixed temperature T. At the left panels of Fig. [3] we 
plot how QD and EoF depend on the anisotropy of the 
J x and J y coupling constants for J z fixed. We note that 
after a certain T the derivatives of QD with respect to 
A = J x — J y are undefined at A = 0. This does not 
happen to EoF. The discontinuity of the derivative shown 
here is similar to what was observed to the dynamical 
decay rate of QD under decoherence in Ref. [12], where 
the authors coined the name "sudden change" for this 
behavior [22] . At the right panels of Fig. [3] we see the 
sudden change of QD keeping all parameters fixed but J z . 
As before, we do not see sudden change for the EoF. The 
sudden change of QD with the coupling constants here 
and the one with time in [T2] suggest that this behavior 
may be related to a technicality of QD rather than to a 
change in a physical property of the system. However, 
no proof of this point is given and further investigation 
is needed. 
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FIG. 3: (Color online) QD (a) and EoF (b) as a function of 
A = J x — J y for J x + J y = 2, J z = 1, and various values of 
kT. QD (c) and EoF (d) as a function of J z for J x + J y = 1, 
A = 7, and various values of kT. For the dotted/black line 
kT = 0.01, dashed/red line kT = 0.1, dash-dotted/blue line 
kT = 0.6, and solid/orange line kT = 1. 

When the magnetic field B is not zero we first analyze 
what happens to the Ising model (J x = J and J y = 
J z =0). Looking at Fig. [I] we see that for the regions 
where EoF is zero QD is negligible. However, as well as 
with EoF [9], QD initially increases as we increase the 
value of 5, going to zero with increasing field. This is 
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true since the density operator Q is |11) (11| (separable 
state) when B — >> oo. On the other hand, as we will 
shortly see, this is not a general result. The behavior of 
QD and EoF can be quite different from each other if we 
work with another model. 




FIG. 4: (Color online) EoF (a) and QD (b) as a function of 
the magnetic field B and the absolute temperature kT for the 
Ising model when J = 1. 

For the XY model (J x , J y ^ and J z = 0) with a 
transverse magnetic field (5^0), for example, QD be- 
haves in quite different (and interesting) ways from what 
we have seen for the Ising model. Also, for the XY model 
we see a new effect, the regrowth of QD. Contrary to the 
behavior of EoF, where we see its sudden death and then 
a revival [9], for QD there is no sudden death. See Fig. 
[5] QD decreases with T, retaining appreciable values and 
then after a critical temperature T c it starts increasing 
again. This is what we call regrowth. Note that in the 
behavior of EoF we never see a regrowth. Indeed, EoF 
increases after decreasing with T only after reaching the 
value zero (sudden death). 









(b) 










(d) 




kT 


kT 



FIG. 5: (Color online) QD, panels (a) and (b), and EoF, 
panels (c) and (d), as a function of the absolute temperature 
kT for the XY model with transverse magnetic field B. For 
(a) and (c) J x — J y — 1 and for (b) and (d) J x = 1.3 and J y = 
0.7. The values for B are 1.1 (red/solid line), 2.0 (blue/dashed 
line), and 2.5 (black/dotted line). 

Moreover, looking at (a) and (c) of Fig. [5] we see 
that QD becomes non-null before the appearance of EoF. 
Also, QD continues to be non-null after EoF disappears. 
If we now analyze (b) and (d) of Fig. [5] we notice that we 
have for all three curves regimes in which QD increases 
while EoF decreases. This is a quite remarkable behavior 
since the decrease of a certain quantum aspect (entangle- 



ment) is simultaneous to the increase of another quan- 
tum aspect (quantum correlations), illustrating clearly 
the distinctive aspects of these two concepts. 

Conclusions. We examined the behavior of the quan- 
tum discord (QD) for a pair of qubits described by the 
Heisenberg model in thermal equilibrium with a reservoir 
at temperature T. By changing the temperature and also 
by applying an external magnetic field we observed sev- 
eral remarkable effects for QD, many of them in sharp 
contrast to the behavior observed for the entanglement 
of formation (EoF) between the two qubits. We found 
that for the XXX model QD signals a quantum phase 
transition for finite T while EoF does not. Also, we ob- 
served regimes where QD increases while EoF decreases 
with T. Moreover, and surprisingly, we showed that for 
the XXZ model there exist regions in the parameter space 
in which EoF is zero while the classical correlation de- 
creases and only QD increases with T. Finally, we also 
observed a new effect for QD which we called regrowth: 
QD decreases with T and starts increasing again after 
reaching a minimum different from zero. 

T.W. thanks CNPq for partially funding this research. 
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